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Salinity  Processes  in  the  Upper-ocean  Regional  Study 


Three-Dimensional  Dynamics 

of  Freshwater  Lenses  in  the 
Ocean’s  Near-Surface  Layer 


By  Alexander  V.  Soloviev, 
Silvia  Matt, 
and  Atsushi  Fujimura 


ABSTRACT.  Convective  rains  in  the  Intertropical  Convergence  Zone  produce 
lenses  of  freshened  water  on  the  ocean  surface.  Due  to  significant  density  differences 
between  the  freshened  and  saltier  seawater,  strong  pressure  gradients  develop,  resulting 
in  lateral  spreading  of  freshwater  lenses  in  the  form  of  gravity  currents.  Gravity 
currents  inherently  involve  three-dimensional  dynamics.  As  a  type  of  organized 
structure,  gravity  currents  may  also  interact  with,  and  be  shaped  by,  the  ambient 
oceanic  and  atmospheric  environment.  Among  the  important  environmental  factors 
are  background  stratification  and  wind  stress.  Under  certain  conditions,  a  resonant 
interaction  between  a  propagating  freshwater  lens  and  internal  waves  in  the  underlying 
halocline  (the  barrier  layer)  may  develop,  while  interaction  with  the  wind  stress  may 
produce  an  asymmetry  in  the  freshwater  lens  and  associated  mixing.  These  two  types 
of  interactions  working  in  concert  may  explain  the  series  of  sharp  frontal  interfaces 
observed  in  association  with  freshwater  lenses  during  the  Tropical  Ocean  Global 
Atmosphere  Coupled  Ocean  Atmosphere  Response  Experiment  (TOGA  COARE).  We 
conducted  a  series  of  numerical  simulations  using  computational  fluid  dynamics  tools. 
These  numerical  experiments  were  designed  to  elucidate  the  relationship  between 
vertical  and  horizontal  fluxes  of  salinity  under  various  environmental  conditions  and 
the  potential  impact  of  these  fluxes  on  the  barrier  layer  and  Aquarius  and  Soil  Moisture 
and  Ocean  Salinity  (SMOS)  satellite  image  formations. 
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INTRODUCTION 

Buoyancy- driven  surface  currents, 
such  as  propagating  rain-formed  lenses 
or  plumes  (we  use  these  terms  inter¬ 
changeably),  are  an  important  compo¬ 
nent  of  the  tropical  ocean  environment 
in  the  Intertropical  Convergence  Zone 
(ITCZ;  see,  e.g.,  Wijesekera  et  al.,  1999). 
These  currents  contribute  to  water  mass 
exchange  by  horizontal  advection  and 
enhanced  vertical  mixing.  As  buoyancy- 
driven  flows,  they  are  a  type  of  organized 
structure  that  resembles  a  classical  grav¬ 
ity  current.  The  water  flow  in  the  lead¬ 
ing  edge  of  the  gravity  current  and  trail¬ 
ing  fluid  contains  a  complex  pattern  of 
three-dimensional  motions  (Ozgokmen 
et  al.,  2004).  Gravity  currents  may  also 
interact  with  ambient  stratification  in  a 
resonant  way,  leading  to  fragmentation 
of  the  near- surface  lens,  and  with  the  sur¬ 
face  wind  stress,  resulting  in  lens  asym¬ 
metry  in  structure  and  mixing  (Simpson, 
1987;  Soloviev  and  Lukas,  1997;  Matt 
et  al.,  2014).  The  patterns  in  gravity  cur¬ 
rents  are  inherently  three  dimensional. 

Notably,  the  dynamics  of  the  fresh¬ 
water  lenses  produced  by  convective 
rains  can  be  linked  to  larger- scale  fea¬ 
tures  (e.g.,  barrier  layer  and  salinity 


fronts),  thus  influencing  global- scale  pro¬ 
cesses  and  contributing  to  the  salinity 
field  detected  in  the  Aquarius  and  Soil 
Moisture  and  Ocean  Salinity  (SMOS)  sat¬ 
ellite  footprints.  The  barrier  layer  (Lukas 
and  Lindstrom,  1991)  separates  the  well- 
mixed  surface  layer  from  the  thermo  - 
dine,  which  enhances  the  air- sea  coupled 
response.  A  substantial  amount  of  data  on 
freshwater  plumes  was  collected  during 
the  Tropical  Ocean  Global  Atmosphere 
Coupled  Ocean  Atmosphere  Response 
Experiment  (TOGA  COARE).  Soloviev 
and  Lukas  (2014)  summarize  these  mate¬ 
rials.  The  NASA  Salinity  Processes  in  the 
Upper-ocean  Regional  Study  (SPURS) 
obtained  new  and  interesting  data  on 
freshwater  plumes  (e.g.,  Anderson  and 
Riser,  2014;  Asher  et  al.,  2014).  However, 
these  data  are  still  fragmentary;  more 
details  about  freshwater  plume  dynamics 
are  required  before  these  subgrid- scale 
phenomena  can  be  incorporated  into 
global  ocean  circulation  models  and 
operational  algorithms  for  sea  surface 
salinity  satellites. 

Oceanic  advection  and  mixing  funda¬ 
mentally  affect  the  sea  surface  salinity  sig¬ 
nal  sensed  by  the  two  satellites.  Improved 
understanding  of  these  processes  is  one 


of  the  SPURS  goals  (Farrar  et  al.,  2014). 
In  this  paper,  we  apply  computational 
fluid  dynamics  tools  to  investigate  hor¬ 
izontal  spreading  and  vertical  mixing 
of  the  freshwater  plumes  produced  by 
convective  rains  in  the  ITCZ.  We  then 
briefly  explain  the  model  setup  and 
describe  the  computational  experiments 
designed  to  illustrate  the  interaction  of 
rain-formed  plumes  with  their  environ¬ 
ments.  The  modeling  results  are  com¬ 
pared  with  conceptual  models  developed 
during  TOGA  COARE. 

MODEL  DESCRIPTION 

Because  the  processes  involved  in  the 
propagation  and  mixing  associated  with 
near-surface  freshwater  plumes  are  three 
dimensional,  explicitly  modeling  their 
dynamics  requires  a  high-resolution, 
nonhydro  static  three-dimensional  mod¬ 
eling  approach.  Thus,  to  simulate  the 
low- density  plume,  we  implemented  a 
computational  fluid  dynamics  model 
with  the  numerical  domain  shown  in 
Figure  la,  following  Matt  et  al.  (2014). 
Periodic  boundary  conditions  are  set  on 
the  walls  along  the  tank  direction.  The 
boundary  conditions  at  the  bottom  and 
the  side  walls  of  the  numerical  tank  are 
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FIGURE  1.  Numerical  domain 
and  model  initialization  (red  color 
corresponds  to  lower  salinity). 

(a)  500  m  x  40  m  x  40  m  domain 
containing  3,120,000  cells.  Az  =  1  cm 
(at  the  top).  Ax  =  50  cm,  Ay  =  1  m. 

(b)  A  close-up  of  the  initial  condi¬ 
tions  showing  the  low  salinity  anom¬ 
aly  produced  by  convective  rain 
at  the  top  of  the  domain,  with  the 
halocline  below. 
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specified  as  zero  shear,  and  zero  heat  and 
salinity  fluxes  are  set  at  the  bottom,  side 
walls,  and  top  of  the  tank.  We  initial¬ 
ize  the  model  with  a  rain  cell  (Figure  la) 
simulating  randomly  distributed  rainfall 
near  the  top  of  the  domain  (Figure  lb). 
Salinity  stratification,  for  example,  a  halo- 
cline  representing  the  barrier  layer  (Lukas 
and  Lindstrom,  1991),  can  be  added  to 
the  model  during  initialization,  and  wind 
stress  can  be  applied  at  the  top  of  the 
tank.  We  tested  a  number  of  wind  stress 
and  halocline  parameter  combinations. 


In  this  work,  we  used  a  large  eddy  sim¬ 
ulation  (LES)  to  model  freshwater  lens 
dynamics.  Two  LES  algorithms  were 
implemented:  LES  with  Wall- Adapting 
Local  Eddy-Viscosity  (WALE)  as  a 
subgrid-scale  model  (Nicoud  and  Ducros, 
1999)  and  a  Hybrid  LES-Detached  Eddy 
Simulation  (DES)  (Strelets,  2001).  Both 
models  have  advantages  and  disadvan¬ 
tages.  WALE  is  optimized  for  rigid-lid 
boundary  conditions  (e.g.,  bottom  and 
side  walls  in  our  simulations);  however, 
it  results  in  excessive  surface  drift  veloc¬ 
ity  when  wind  stress  is  applied  at  the  top 
of  the  numerical  tank  due  to  absence  of 
surface  waves.  DES  results  in  more  real¬ 
istic  surface  drift  (Matt  et  al.,  2011)  but 
produces  less  realistic  results  near  the 
rigid  walls  where  the  specified  zero- 
shear  boundary  condition  is  applied. 


In  this  article,  depending  on  the  case, 
we  use  either  the  LES  WALE  or  DES 
model.  More  details  on  the  numerical 
schemes  and  setup,  including  the  WALE 
and  DES  formulations,  can  be  found  in 
Matt  et  al.  (2011,2014). 

NUMERICAL  SIMULATIONS 
AND  RESULTS 

We  conducted  high-resolution  numer¬ 
ical  experiments  on  the  dynamics  of 
near- surface  buoyant  plumes  in  the  pres¬ 
ence  of  ambient  stratification  and  wind 


stress.  The  following  process  studies  were 
performed  using  the  computational  fluid 
dynamics  model  described  in  the  pre¬ 
vious  section.  Depending  on  the  case, 
environmental  conditions  for  the  evolv¬ 
ing  buoyant  plume  are  set  either  via  ini¬ 
tial  conditions  (salinity  stratification)  or 
as  a  boundary  condition  at  the  top  of  the 
domain  (wind  stress). 

Freshwater  Lens  Propagating 
Under  Calm  Weather  Conditions 
in  Nonstratified  and  Stratified 
Environments 

Figure  2  shows  the  spreading  of  a  rain- 
formed  plume  in  an  unstratified  envi¬ 
ronment,  simulated  using  the  LES  WALE 
algorithm.  The  low- salinity  lens  was  ini¬ 
tialized  in  the  near- surface  layer  as  a  salin¬ 
ity  anomaly  within  a  5  m  deep  box  with 


horizontal  dimensions  of  40  m  x  40  m  in 
the  x  and  y  directions.  The  initial  salin¬ 
ity  anomaly  was  randomly  distributed  at 
a  1.2  psu  peak  and  0.6  psu  average  level, 
well  within  the  range  of  near- surface 
salinity  anomalies  observed  in  the  west¬ 
ern  equatorial  Pacific  (Soloviev  and 
Lukas,  2014). 

As  the  plume  evolved  from  its  ini¬ 
tial  state  (see  Figure  1),  randomly  dis¬ 
tributed  rain  signatures  quickly  merged 
into  a  continuous  low- salinity  plume. 
Figure  2a  shows  contour  plots  of  the 
freshwater  plumes  velocity  field  at  the 
surface  and  at  13  m  depth  after  2,500  s  of 
simulation.  The  simulation  was  limited 
to  2,500  s  because  the  spreading  plume 
was  approaching  the  domains  boundar¬ 
ies.  The  frontal  edges  of  the  plume  spread 
with  a  speed  of  0.14  m  s_1  along  the  tank. 
The  modeled  current  velocities  are  of  the 
same  order  of  magnitude  as  observations 
(Wijesekera  et  al.,  1999;  Soloviev  et  al., 
2002).  The  spreading  plume  induces 
velocity  fluctuations  of  the  order  of 
0.004  m  s-1  below  the  plume  fronts  at  a 
13  m  depth.  Because  the  initial  density 
anomaly  was  in  the  center  of  the  tank, 
the  velocity  pattern  in  this  simulation 
was  almost  symmetrical  with  respect  to 
the  center  of  the  tank.  After  2,500  s  of 
simulation,  the  freshwater  lens  signifi¬ 
cantly  expanded  horizontally  but  col¬ 
lapsed  in  the  vertical  direction,  from  5  m 
to  roughly  1  m  thickness  (Figure  2d).  It 
is,  however,  necessary  to  remember  that 
this  simulation  was  performed  for  calm 
weather  conditions. 

The  contour  plot  of  the  salinity  gra¬ 
dient  in  Figure  2b  reveals  a  streak-like 
structure  in  the  fluid  trailing  behind  the 
plume  edge.  The  streaks  are  approxi¬ 
mately  aligned  with  the  plume  propaga¬ 
tion  direction  but  are  intermittent.  The 
cross-section  contour  plot  of  the  verti¬ 
cal  velocity  shown  in  Figure  2c  suggests 
that  these  structures  have  the  form  of 
rolls.  The  nature  of  these  streaks  in  the 
ocean  is  still  not  completely  clear;  they 
may  be  of  the  same  nature  as  Lesieur  s 
(2008)  streaks  or  Kuettners  (1971)  atmo¬ 
spheric  rolls.  These  coherent  structures 


These  numerical  experiments  were  designed 
to  elucidate  the  relationship  between  vertical 
and  horizontal  fluxes  of  salinity  under  various 
environmental  conditions  and  the  potential  impact 
of  these  fluxes  on  the  barrier  layer  and  Aquarius 
and  Soil  Moisture  and  Ocean  Salinity  (SMOS) 
satellite  image  formations. 
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can  be  caused  by  vertical  shear  within  the 
propagating  plume.  Note  that  Ozgokmen 
et  al.  (2004)  reproduced  similar  coherent 
structures  in  a  three-dimensional  numer¬ 
ical  simulation  of  an  overflow  gravity  cur¬ 
rent.  The  width  of  streaks  is  an  order  of 
magnitude  larger  than  the  thickness  of 
the  freshwater  lens,  and  their  length  is 
two  orders  of  magnitude  larger  than  the 
thickness  of  the  freshwater  lens.  These 
numbers  are  consistent  with  Lesieurs 
(2008)  conceptual  model. 

With  the  addition  of  the  halocline 
below  the  initial  freshwater  anomaly,  fur¬ 
ther  interaction  was  observed.  In  this 
experiment,  we  initialized  the  halocline 
with  the  linear  salinity  change  of  4  psu 
between  6  m  and  16  m  water  depth.  This 
initialization  resembles  a  halocline  (bar¬ 
rier  layer)  reported  by  Vinayachandran 
et  al.  (2002)  in  the  northern  part  of  the 
Bay  of  Bengal. 

Figure  3  shows  contour  plots  of  the 
velocity  field  on  the  surface  and  at  13  m 
depth  after  2,500  s  of  simulation.  In  the 
case  of  no  stratification,  the  propagating 
freshwater  plume  induced  the  circulation 
below  its  edge  (Figure  2a).  In  the  case  of 
environmental  stratification  (halocline), 
plume  propagation  resulted  in  the  exci¬ 
tation  of  internal  waves  in  the  underlying 
halocline  and  apparent  modulation  of 
the  lens  structure  by  these  internal  waves 
(Figure  3a).  The  internal  waves  excited  in 
the  barrier  layer  by  the  freshwater  lens 
propagate  somewhat  ahead  of  the  lens 
edge  (Figure  3c). 

Based  on  the  laboratory  results  of 
Simpson  (1987),  Soloviev  and  Lukas 
(1997)  hypothesized  that  the  freshwater 
lens  in  the  oceans  near-surface  layer  may 
interact  with  environmental  stratifica¬ 
tion.  A  small  thickness  of  the  near-surface 
anomalies  is  favorable  for  the  develop¬ 
ment  of  nonlinear  interactions.  These 
laboratory  studies  indicate  that  any  den¬ 
sity  anomaly  in  the  upper  layer  tends  to 
spread  horizontally.  The  leading  edge 
of  the  gravity  current  along  a  boundary 
through  a  uniformly  stratified  medium 
can  generate  several  different  modes 
of  internal  waves. 
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FIGURE  2.  (a)  The  velocity  magnitude  field  at  the  surface  and  at  13  m  depth  induced  by  the  fresh¬ 
water  lens  spreading  from  the  initial  position  of  the  plume  at  the  center  of  the  domain  into  a 
non-stratified  environment  after  2,500  s  of  simulation,  (b)  Coherent  structures  seen  in  the  magni¬ 
tude  of  the  salinity  gradient  on  the  top  of  the  spreading  freshwater  plume,  (c)  Contour  plot  of  the 
vertical  velocity  component  reveals  rolls  inside  the  propagating  plume,  (d)  Salinity  and  magnitude 
velocity  contour  plots  near  the  plume  edge  (side  view). 
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According  to  the  laboratory  observa¬ 
tion  described  by  Simpson  (1987),  “These 
waves  affected  the  form  of  the  gravity 
current  behind  the  head  in  a  rhythmi¬ 
cal  manner.  The  fluid  in  the  original  head 
was  cut  off  from  the  flow,  which  formed 
a  second  head.  The  process  was  repeated 
and  later  a  third  new  front  appeared...” 
However,  this  striking  effect  was 
observed  only  when  the  following  con¬ 
ditions  for  the  Froude  number  (Fr)  and 
the  fractional  depth  of  the  gravity  current 
(r)  were  satisfied: 


and 

r=  h/H  <0.2  (2) 

where  U  is  the  speed  of  the  gravity  cur¬ 
rent,  h  is  the  thickness  of  the  gravity  cur¬ 
rent,  H  is  the  depth  of  the  ocean,  N  is  the 
Brunt- Vaisala  frequency  of  the  strati¬ 
fied  surroundings,  and  n  =  3.14.  This  is 
a  resonant  type  mechanism  of  interaction 
between  the  freshwater  lens  and  the  strat¬ 
ified  environment. 

In  the  case  shown  in  Figure  3,  Fr  =  0.08 
and  r  =  0.05;  according  to  conditions  (1) 
and  (2),  resonant  interaction  between 
the  freshwater  plume  and  the  under¬ 
lying  halocline  was  expected.  For  the 
no -stratification  case  shown  in  Figure  2, 
there  is  no  halocline  (Fr  =  oo),  so  no  reso¬ 
nant  interaction  was  possible. 

The  banding  pattern  on  the  surface 
(Figure  3a)  can  be  interpreted  as  an  indi¬ 
cation  of  the  plume  fragmentation  due 
to  resonant  interaction  with  the  ambient 
stratification.  Notably,  this  fragmenta¬ 
tion  interrupted  the  coherent  rolls  as  well 
(Figure  3b). 


FIGURE  3.  (a)  The  velocity  magnitude  field  at  the  surface  and  at  13  m  depth  induced  by 
the  freshwater  lens  propagating  above,  and  interacting  with,  the  underlying  halocline. 
(b)  Fragmentation  of  the  plume  into  bands  due  to  resonant  interaction  with  the  underlying 
halocline.  (c)  Contour  plots  of  the  velocity  magnitude  (side  view)  showing  generation  of  inter¬ 
nal  waves  by  the  spreading  lens.  Note  the  internal  waves  propagating  ahead  of  the  lens  edge. 


Interaction  of  Freshwater  Plume 
with  Wind  Stress 

Next,  we  investigated  the  case  of  a  fresh¬ 
water  plume  interacting  with  wind 
stress  in  a  nonstratified  environment. 
In  order  to  produce  a  realistic  aque¬ 
ous  viscous  sublayer,  a  DES  scheme  was 
implemented.  We  applied  wind  stress 
at  the  top  of  the  numerical  tank  and 
spun  up  the  model  for  500  s.  The  plume 
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simulation  was  then  initialized  with  a 
low- salinity  plume,  imitating  randomly 
distributed  rainfall  similar  to  Figure  1 
but  with  no  stratification  below  the  fresh¬ 
water  plume.  We  tested  the  three  cases 
corresponding  to  wind  speeds  at  10  m 
height,  Ui0  -  0  m  s-1,  U10  =  4.0  m  s-1,  and 
U10  =  8.5  m  s_1  (Figure  4).  The  plume 
structure  showed  a  notable  asymmetry 
relative  to  the  wind  direction  when 
forced  with  wind  stress.  The  asym¬ 
metry  of  the  plume  apparently  increases 
with  wind  speed. 

The  plume  asymmetry  caused  by  wind 
stress  is  clearly  seen  in  Figure  5,  which 
shows  two  isosurfaces  associated  with  the 
freshwater  plume.  Figure  6  shows  evo¬ 
lution  of  the  freshwater  plume  in  time. 
Notably,  mixing  on  the  upwind  side  of  the 
plume  is  more  intense  than  on  the  down¬ 
wind  side  of  the  plume. 

The  asymmetry  of  freshwater  plumes 
relative  to  the  wind  speed  direction 
was  observed  in  the  western  equatorial 
Pacific  during  TOGA-COARE  (Soloviev 
and  Lukas,  1997;  Soloviev  et  al.,  2002). 
Soloviev  and  Lukas  (1997)  suggested 
that  this  asymmetry  can  be  explained 
by  the  mechanism  of  Stommels  (1993) 
Overturning  Gate,  schematically  shown 
in  Figure  7.  When  the  wind  stress  is 
directed  toward  lower  density,  gravi¬ 
tational  instability  may  trigger  Kelvin- 
Helmholtz  instability  and  intensive  ver¬ 
tical  and  cross -frontal  mixing  (Soloviev 
and  Lukas,  1997).  In  the  process  of  mix¬ 
ing,  the  interface  may  either  reach  the 
compensated  state  in  density  (when 
the  temperature  differences  associated 
with  the  low- salinity  lens  are  accounted 
for)  or  entirely  disappear;  in  both  cases, 
the  interface  is  no  longer  prominent  in 
the  density  field.  When  the  wind  stress 
is  directed  toward  higher  density,  sta¬ 
ble  stratification  increases  due  to  lat¬ 
eral  advection  of  less  dense  water,  which 
inhibits  mixing.  Consequently,  the 
downwind  frontal  interface  “freezes”  and 
can  drift  downwind  until  the  wind  sub¬ 
stantially  changes  either  in  direction  or 
speed  with  respect  to  the  plume. 

According  to  the  criteria  proposed 


by  Soloviev  et  al.  (2002),  instability  on 
the  upwind  edge  of  the  plume  can  take 
place  when 

ReL~y601l3(g’L/uiy13  <2.74  (3) 

where  g  =  gAp/p  is  the  reduced  gravity, 
Ap  the  density  difference  across  the  lead¬ 
ing  edge  of  the  plume,  p  the  water  den¬ 
sity,  w*  the  frictional  velocity  in  water,  and 
the  L  leading  edge  width  (the  distance  to 
which  the  denser  water  propagates  over 
the  plume  edge).  For  the  case  shown  in 
Figure  6,  ReL  ~  1,  which  satisfies  inequal¬ 
ity  (3).  There  are  indeed  indications  of 
intense  mixing  on  the  upwind  side  of  the 
plume  (Figure  6),  and  this  is  consistent 


with  theory.  Note  that  freshwater  plumes 
with  a  larger  density  difference  and/or 
lower  wind  speed  may  not  satisfy  condi¬ 
tion  (3).  This  theory  also  ignores  Coriolis 
forces  and,  strictly  speaking,  is  valid  only 
for  the  equatorial  region. 

CONCLUSIONS 

Convective  rains  within  the  ITCZ  pro¬ 
duce  localized  freshwater  plumes,  which 
tend  to  spread  horizontally.  These  plumes 
resemble  gravity  currents  and  can  interact 
with  background  stratification  and  wind 
stress.  The  dynamics  of  these  near-surface 
plumes  are  inherently  three  dimensional, 
and  explicit  modeling  of  the  processes 
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FIGURE  4.  Contours  of  the  near-surface  freshwater  plume  interacting  with  wind  stress  (side  view). 
The  length  scale  is  the  same  for  both  coordinates. 
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FIGURE  5.  A  bird’s-eye  view  of  the  freshwater  lens  under  wind  stress  action  after  a  1,000  s  simula¬ 
tion  contoured  with  salinity  isosurfaces  at  34.92  psu  (brown)  and  34.97  psu  (pink).  The  length  scale 
is  the  same  for  all  three  coordinates. 
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FIGURE  6.  Evolution  of  the  freshwater  plume  interacting  with  wind  stress.  The  contour  plot  of  salinity  (side  view)  reveals 
strong  asymmetry  of  the  freshwater  lens  structure  and  mixing  relative  to  wind  direction.  Note  that  the  red  color  here  cor¬ 
responds  to  lower  salinity.  The  length  scale  is  the  same  for  horizontal  and  vertical  axes.  A  related  animation  is  available  in 
the  online  supplementary  material  at  http://www.tos.org/oceanography/archive/28-1_soloviev.html. 
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FIGURE  7.  Depiction  of  Stommel’s  Overturning  Gate  in  application  to  the  near-surface  freshwater 
plume  propagating  as  the  gravity  current  and  interacting  with  wind  stress,  as  proposed  by  Soloviev 
and  Lukas  (1997). 


involved  requires  high-resolution,  non¬ 
hydrostatic  three-dimensional  models. 
In  this  work,  we  use  computational  fluid 
dynamics  to  reproduce  the  interaction  of 
a  freshwater  lens  with  a  stratified  envi¬ 
ronment  and  wind  stress.  Our  model 
reproduces  development  of  coherent 
structures  (streaks  or  rolls)  in  the  prop¬ 
agating  plume  and  resonant  interaction 
between  a  near- surface  freshwater  plume 
and  ambient  stratification,  which  can  lead 
to  repeating  frontal  interfaces  or  “surface 
lines”  that  have  been  previously  reported 
from  field  observations  (Soloviev  and 
Lukas,  1997).  However,  coherent  streaks 
await  detailed  investigation. 

Our  results  also  indicate  that  wind 
stress  interacting  with  a  near-surface 
freshwater  plume  can  lead  to  marked 
asymmetry  in  sea  surface  signature,  prop¬ 
agation,  and  mixing  associated  with  the 
plume.  This  is  consistent  with  substantial 


statistics  collected  during  TOGA  COARE 
(Soloviev  et  al.,  2002). 

One  interesting  result  that  can  be  use¬ 
ful  for  planning  future  fieldwork  is  that 
freshwater  lenses  spread  rapidly  horizon¬ 
tally  and  at  the  same  time  collapse  ver¬ 
tically.  This  process  results  in  the  pres¬ 
ence  of  shallow,  low- salinity  layers  in  the 


upper  few  meters  of  the  ocean  under  low 
and  moderate  wind  speed  conditions. 
Soloviev  and  Lukas  (1997)  observed 
such  shallow  freshened  layers  during 
TOGA  COARE.  Apparently,  the  relation¬ 
ship  between  the  vertical  and  horizontal 
structures  of  near- surface  salinity  fields 
depends  on  environmental  conditions 
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(wind  and  waves),  and  a  full  investigation 
is  planned  in  the  Pacific  during  SPURS-2. 
These  processes  potentially  impact  bar¬ 
rier  layer  physics  and  thus  affect  Aquarius 
and  SMOS  satellite  image  formation. 

Influx  of  freshwater  due  to  convec¬ 
tive  rains  can  result  in  substantial  local¬ 
ized  salinity  anomalies  (freshwater  lenses) 
in  the  upper  few  meters  of  the  ocean. 
The  Aquarius  and  SMOS  missions  do 
not  explicitly  resolve  fine  surface  salinity 
structures  because  the  spatial  resolution 
of  these  satellites  is  of  the  order  of  100  km. 
Nevertheless,  understanding  the  mixing 
and  advection  processes  associated  with 
freshwater  lenses  is  important  for  calibra¬ 
tion  and  correct  interpretation  of  satellite 
salinity  imagery  (Farrar  et  al.,  2014). 

Future  modeling  work  should  investi¬ 
gate  the  dynamics  of  freshwater  plumes 
under  the  influence  of  both  ambient 
stratification  and  wind  stress  and  how 
they  interact  to  affect  plume  dynam¬ 
ics.  The  Coriolis  force  and  tempera¬ 
ture  anomalies  in  the  freshwater  plume 
should  be  included  in  future  models  as 
well.  Synthetic  aperture  radar  data,  which 
under  certain  conditions  can  resolve  hor¬ 
izontal  features  on  the  sea  surface  associ¬ 
ated  with  freshwater  plumes  down  to  the 
1  m  scale,  should  be  useful  for  model  ver¬ 
ification  (Matt  et  al.,  2014).  M 

ONLINE  SUPPLEMENT.  A  related  animation  of 
the  evolution  of  a  freshwater  plume  interacting  with 
wind  stress  is  available  in  the  online  supplemen¬ 
tary  material  at  http://www.tos.org/oceanography/ 
archive/28-1_soloviev.html. 
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